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Differentially rotating disks of dust 

Marcus Ansor^ and Reinhard MeinelQ 



Abstract 

We present a three-parameter family of solutions to the stationary 

axisymmetric Einstein equations that describe differentially rotat- 

Q ■ ing disks of dust. They have been constructed by generalizing the 



Neugebauer-Meinel solution of the problem of a rigidly rotating disk 
of dust. The solutions correspond to disks with angular velocities de- 
pending monotonically on the radial coordinate; both decreasing and 



xjTi ' increasing behaviour is exhibited. 

^ ■ In general, the solutions are related mathematically to Jacobi's in- 

(^ . version problem and can be expressed in terms of Riemann theta 
functions. A particularly interesting two-parameter subfamily repre- 

Qs^ I sents Backlund transformations to appropriate seed solutions of the 

"^ ' Weyl class. 
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1 INTRODUCTION 



Although many rigorous solutions to Einstein's field equations have been con- 
structed, only few of them can be applied to physically relevant situations. In 
the particular field of stationary axisymmetric solutions so far only the Kerr black 
holes |]l| and the rigidly rotating disks of dust @, ^ §] are known to be physically 
relevant. On the other hand, stationary axisymmetric solutions are of special in- 
terest to astrophysics to model equilibrium configurations like stars and galaxies. 
The complicated structure of the Einstein equations describing the interior of a 
rotating body gives httle hope for a rigorous global solution in the near future. 
However, if one restricts oneself to considering only dust configurations of finite 
extension, the body fiattens in an extreme manner (it becomes a disk) and the 
interior equations turn into boundary conditions for the exterior vacuum equa- 
tions. Now these equations can be expressed by a single nonlinear equation - 
the so-called Ernst equation P, ^. For treating the Ernst equation, there are 
powerful analytic methods available which come from soliton theory [|^-|12[. In 



particular, Korotkin |T^, Q and Meinel and Neugebauer |jT5|, see also [|l^, |T 



were able to construct a class of solutions containing a finite number of complex 
parameters and one arbitrary real solution to the axisymmetric three dimensional 
Laplace equation. In this paper we present a three-parameter subclass describing 
disks of dust revolving with a non uniform angular velocity. These solutions are 
analytic in the sense that they belong to the class of solutions under discussion 
and therefore strictly satisfy the Ernst equation. On the other hand, these solu- 
tions are numerical solutions since the real-valued potential function mentioned 
above has been determined numerically in order to satisfy the boundary condi- 
tions. The accuracy that has been obtained was very high (generally 12 digits) 
such that for any practical use the solutions are just as good as purely analytic 
ones. Moreover, the accuracy may in principle be increased arbitrarily. 
The paper is organized as follows. In the first chapter, the boundary value prob- 
lem for differentially rotating disks of dust is introduced and the class of solutions 
in question is reviewed. The numerical methods by which we were able to obtain 
our subclass of differentially rotating disks of dust will be discussed in the first 
part of the second chapter. This will be followed by a thorough illustration of the 
parameter space of these solutions. The subsequent subchapters contain detailed 
discussions about particular limits. 

In what follows, units are used in which the velocity of light as well as Newton's 
constant of gravitation are equal to 1. 



1.1 Metric Tensor, Ernst equation, and boundary condi- 
tions 

The metric tensor for axisymmetric stationary and asymptotically flat space- 
times reads as follows in Weyl-Papapetrou-coordinates (p, C,V^, t): 

ds" = e-2^[e2^(dp2 + d^) + p'd^''] - e^^'idt + a d^f . 

For this line element, the vacuum field equations are equivalent to a single com- 
plex equation - the so-called Ernst equation^ 

(5R/) A/ = (V/)^ (1) 

A-— -— — v-(— — 

dp^ pdp 9C^' \dp^ d( 



where the Ernst potential / is given by 

- — a,p, b^p = 

P ^ P 



f = e^^ + ib with &,(^ = a,p, &,p = a,c- (2) 



To obtain the boundary conditions for differentially rotating disks of dust, one 
has to consider the field equations for an energy-momentum-tensor 

T'' = eu'u^ = ap{p)e^-''6{C)u'u'', 

where e and dp stand for the energy- density and the invariant (proper) surface 
mass-density, respectively, 6 is the usual Dirac delta-distribution, and m* denotes 
the four- velocity of the dust material^. 

Integration of the corresponding field equations from the lower to the upper side 
of the disk (with coordinate radius po) yields for ^ = 0+ and < p < po the 
conditions (see |jl8|) 



2ncx, = e^'-^iU^ + ^Q) (3) 

e^V + Q(e'^),C + (M' = (4) 

with 

g = -pe-^^[6>,C + (e2^),,(e^^),d. (5) 

Note that boundary condition (^ for the Ernst potential / does not involve the 
surface mass-density Up. This condition comes from the nature of the material the 



^ The remaining function k can be calculated from the Ernst potential / by a line integration, 
■^u' has only (p- and t- components. 



disk is made of. Therefore, equation (^) will be referred to as the dust- condition. 
The angular velocity Q = u'^/u^ of the disk can be calculated from 

fi = ^^- (6) 

The following requirements resulting from symmetry conditions and asymptotical 
flatness complete our set of boundary conditions: 

• Regularity at the rotation axis is guaranteed by 

|(o.C)^o. 

• At infinity asymptotical flatness is realized by f/ — > and a —>■ 0. For 
the potential b this has the consequence h —^ h^o = const. Without loss of 
generality, this constant can be set to 0, i.e. / ^ 1 at infinity. 

• Finally, we assume refiectional symmetry with respect to the plane C = 0^ 
i.e. /(p, — C) = /(p, C) (with a bar denoting complex conjugation). 

1.2 Solutions related to Jacobi's inversion problem 

Meinel and Neugebauer [|T^] showed that for an arbitrary integer p the function 
/ defined by 




^'^^ KPdK 



W{K) 



Vp 



W{K) 



\ 



{K + iz){K - iz) il{K~ K,){K - K,) 



v=l 



satisfies the Ernst equation^. Hereby, the Ky are arbitrary complex parameters. 
The variable z = p + iC, is the complex combination of the coordinates p and 
C,. The (^-dependent) values for the K'^'^^ as well as the integration paths on a 
two-sheeted Riemann surface have to be taken from the solution to the following 
Jacobian inversion problem: 

" ^"^' K^dK 



v=l 



E/ m^ = ^^' ^<3<P- 



K^ 



W{K) 



^Korotkin [jTsI, pi, see also fl^, 11^, found solutfons to the Ernst equatfon which are cfosely 
related to the solutions considered here. 



The potential functions Vj (0 < j < p) may be any real solutions to the axisym- 
metric Laplace equation /\vj = satisfying the recursion conditions 

i^i,^ = 2^i-i + ^'"j~i,z- (7) 

These recursion conditions are automatically satisfied by the ansatz 
1 f K^H{K) 



Vi = / ' ' ^ dK, (8) 

where S is some curve (or even some set of curves) in the complex plane. With (^, 
the regularity of the resulting solution at (p, Q = {\'^[K,^]\, ^[K,^]) is guaranteed. 
However, there are discontinuities along the curve S' = {(p, C) • C i ip ^ ^j- 
Moreover, a free function H defined on S enters the class of solutions. Hence, 
this ansatz allows us to consider a restricted class of boundary value problems^ 
in which the curve S results from the shape of the boundary S'. 
Ernst potentials with reflectional symmetry /(p, —() = f{p, Q are characterized 
by the following properties: 

• For each parameter K^, there is another parameter i^^ with K^j^ = —Ky. 



• H{±K) = H{K) 

For differentially rotating disks, we can set E = {K : K = ipox, ~1 < x < 1}. 
Thus we get 



I" {jxypx^ ^^^ ^i. = v^(i^^c7p^F+Tp7p^ mzD)<o). 



_ f'iixyhix') 

"j — Po 



In this expression we require the real-valued function h to be analytic on the 
interval [0, 1]. This is necessary for an analytic behaviour of the angular velocity 
fl for all p G [0,po]. 
The additional requirement 

Ml) = (9) 

leads to a surface mass-density Cp of the form 



o"p(p) = i'{p)\j pI - P^ (with ip analytic on [0, po]) 



and therefore ensures vanishing ap at the rim of the disk. The resulting Ernst 
potential / depends on the normalized coordinates (p/pcC/Po)) on the parame- 
ters Xy = K^/pq and functionally on h. 

One obtains the special case of rigid rotation |^ when the following equations 
are all satisfied: 



^The restriction comes from the requirement that the potential functions Vj be real. 



p 



^2 - 1 



Xf = -1 + -, X2 = -Xi 



2, arsinh[/i(l — x^)] 



• h{x') = fi , ^^' ' ^. (10) 

7rv/l + /i2(l-x2)2 

The parameter /i is related to the angular velocity 

and runs on the interval (0, /iq) with fiQ = 4.62966184 . . . For fi <^ 1 one obtains 
the Newtonian limit of the Maclaurin disk. On the other hand, fi —>■ fiQ and 
Po ~* yields the ultrarelativistic limit of the extreme Kerr black hole. 

In this article we explore the subclass p = 2 of the solutions introduced above. It 
will be shown how, for a given value of the complex parameter Xi {Ki = poXi), 
the freedom of the choice of the function h with the property (P) has been used 
to find a solution satisfying the dust-condition (^). It turns out that for each Xi 
within a certain region (a more precise description follows) there is a function h 
such that the resulting Ernst potential can be interpreted as having been created 
by a differentially rotating disk of dust. The accompanying surface mass-density 
and angular velocity may afterwards be calculated according to equations (^ and 

2 DIFFERENTIALLY ROTATING DISKS 

2.1 The numerical scheme 

As already mentioned above, in this paper we consider the class of solutions 
introduced in the previous chapter for the particular case p = 2. Here we can 
prescribe 

1. the coordinate radius po 

2. the complex parameter Xi with 3?(Xi) < and S(Xi) < (without loss of 
generality). Then Ki and K2 follow from Ki = poXi, K2 = —Ki. 

3. a real-valued function h : [0, 1] — > M which is analytic everywhere in [0, 1] 
(i.e. in particular at the boundaries of the interval) and vanishes at the 
upper boundary: h{l) = 0. 
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For such a choice all the requirements stated in chapter 1.1 are satisfied except the 
dust-condition (^Q Now, this condition yields a complicated nonlinear integral 
equation for h: 



< p < po, C = 0+, / = e^^ + 16 = /(p/po, C/Po; X,;h), 



(11) 



Note that the resulting function h depends parametrically on Xi (but not on po). 
With expansions of the functions h and D in Chebyshev-polynomials (this can 
be done since both of them are analytic in [0, 1]) we try to discretize equation 

(0): 

^ 1 

• /i(x^) ^ 2_. hjTj^i{2x'^ ^ 1) ~ 7:hi, Tj{t) = cos[j arccos(r)] 

i=i 

N 

h{l) = => hi = -2^hj 

i=2 



N-l 

• D{x^- Xi; h)^Y. ^i(^i; ht)T,^,{2x^ - 1) - -D^{X,; hu) 
i=i 

In this manner, the nonlinear integral equation (|ll|) is approximated by a finite 
system of nonlinear equations 

Dj{Xi- hk) = il<j<N, 2<k< N). (12) 

The system (0) has been solved numerically by a Newton-Raphson-method. For 
this technique one needs a good initial guess for the solution. Fortunately, the 
values hk are given exactly for the rigidly rotating disks. Therefore, we start 
with an Xi that differs only slightly from those for the rigidly rotating disks [say 
Xf = (—1 + e) + i//i] and take the h^s for the rigidly rotating disks as initial 
values. The newly calculated nearby solution serves then as an initial guess for 
another solution further away from the rigidly rotating disks. Thus we can grad- 
ually explore the whole parameter region of Xi . 

The numerical code written to implement this scheme produces results with ex- 
cellent convergence. For almost all values Xi inside the available parameter 



Additionally one has to ensure that the surface mass-density ap given by equation (pi) is 
positive and finite within [0, po]. Furthermore, the global regularity of the Ernst potential has 
to be checked. Fortunately, our solutions possess these properties. 



region^ the cancellation of the terms in equation (|rTj) up to the 12th digit and 
even beyond has been achieved within the whole range x^ G [0, 1]. The resulting 
Chebyshev-coefficients fall off rapidly (generally A^ = 20 suffices to achieve the 
previously mentioned accuracy of 12 digits) and the resulting function h indeed 
has the desired smooth analytic behaviour. As already mentioned in footnote |^, 
the accompanying surface mass-density turns out to be positive and finite and 
the Ernst potential is regular outside the disk. 

2.2 The parameter region of the solutions and examples 

In the following we discuss the solutions obtained, as a function of the parameter 
X\. For each value outside the hatched region in Figure 1 we have found a 
corresponding solution to the Ernst equation satisfying the dust-condition. What 
follows is a discussion of the details of Figure 1: 

1. Differential rotation: For 3fJ(X^) = —1 we find the rigidly rotating disks. 
On the left hand side of this line (i.e. for 3?(Xf) < —1) the solutions 
turned out to possess an angular velocity f2 which increases with the radial 
coordinate p. On the other side, for 3fJ(Xf ) > — 1, the function fi decreases 
as p increases. 

2. Ergoregions: For values Xf inside the area encompassed by the curve F^;, 
the curve Tu and parts of the curves F^ and Fo-, the corresponding disks 
possess an ergoregion, i.e. a portion of the (p, (^)-space within which the 
function e^^ is negative. 



3. Ultrarelativistic limit: As will be shown in chapter p.3| , any simultaneous 
limit 

• po tends to 

• Xl tends to a value on V^ 

turns out to be an ultrarelativistic limit. In the case of non vanishing values 
for p^ + C^) the resulting / tends to the Ernst potential of an extreme Kerr 
black hole. If, on the contrary, finite values for a/p^ + C^/Po are maintained, 
non asymptotically fiat solutions can be obtained. These results are in 
agreement with a conjecture by Bardeen and Wagoner |19[ . 



4. The Newtonian limit \Xl\ -^ oo.- Here the Ernst potential tends to 1, i.e. 
it describes a Minkowski space. In a given neighbourhood (for large values 
of \Xl\) a post-Minkowskian expansion (with the first coefficient being a 
Newtonian potential) of the Ernst potential can be carried out. One finds 
that the resulting Newtonian coefficient is the gravitational field of a rigidly 
rotating Maclaurin-disk. 



^The exceptions are to be found in narrow stripes along the curve r^ (see next chapter) 



5. The Newtonian limit '^{Xi) -^ 0: For real and positive values of Xf (hence 
Xi = —X2, both real) again the corresponding Ernst potential is equal 
to 1. Here, the Newtonian coefficient of the post-Minkowskian expansion 
describes a disk with decreasing fi(p). 

6. Bdcklund limit: For real and negative values of X^ we get Xi = X2. Then 
the complicated structure of the Ernst potential / simplifies considerably. 
One finds that these solutions can be interpreted as Backlund transforms 
of appropriate seed solutions of the Weyl class. A detailed discussion on 
this subclass is given in chapter 



7. The hatched region and the curve T„: Inside the hatched region no solutions 
have been found satisfying both of the requirements 

• / is regular everywhere outside the disk. 

• The function h is analytic for all x^ G [0, 1]. 

Apart from the ultrarelativistic curve Tu the hatched region is encompassed 
by the curve Fo-. Starting from F$ and ending at F^/, the corresponding 
Ernst potentials describe a transition from Minkowski space to the ultrarel- 
ativistic limit, just as the rigidly rotating disks and the Backlund solutions 
(at Vb) do. Now, all solutions along Fo. possess the property that the 
derivative of the surface mass-density vanishes at the rim of the disk, i.e. 
(Tp{p) = (pq — p^)^/^ i'ip) (with ip analytic in [0,po]) . Surprisingly, this 
physical property coincides with the failure of our numerical method at F^. 
Further investigations are necessary to clarify this coincidence. 

Figure 2 shows representative examples of solutions. From Figure 3 we can get 
an impression of how "relativistic" our solutions are. Here, the physical quantity 
M'^/J is plotted over the available parameter region of Xf with the gravitational 
mass M and the total angular momentum J given by 

M = 2f{Tab-^TgabKedV {T = g^bT'^'), J = - j T^.n'^r^'dV. 

{S is the spacelike hypersurface t = constant with the unit future-pointing normal 
vector n"; the Killingvectors ^" and rj°' correspond to stationarity and axisymme- 
try, respectively.) Note that M and J can also be calculated from the behaviour 
of the Ernst potential at infinity: 



t/ = -— + 0(r-2), 5 = _2 J^ + C(r-3) (p = r sin^, C = r cos^). 



The agreement of the two expressions yields an excellent confirmation of the 
regularity of the solutions. For Figure 3, not only the values M'^/J have been 
determined but also this agreement was checked. 



2.3 Ultrarelativistic limits 

In this chapter we show how ultrarelativistic limits of our solutions can be ob- 
tained if simultaneously 

• Xf tends to a value on F^. 

• the coordinate radius po tends to zero. 

To this end we first consider the limit po —>■ for a fixed value Xf ^ Tu and 
finite values of p^ + C^ 7^ 0. Here, / tends to 1. In a second step, we investigate 
the above simultaneous limit again for p^ + C^ 7^ and find that / tends to 
the Ernst potential of an extreme Kerr black hole. Finally, we indicate how 
further non asymptotically flat solutions can be obtained by fixing finite values 

of ^/^Te/po. 



1. The limit po — ^ for Xf ^ Tu: Similar to the treatment in |2^ we rewrite 



the expression for / and Jacobi's inversion problem in the equivalent form^ 
with 

Vj = Vj - Wj =^J- J yy (j = 0, 1, 2). 

In the limit po ^ 0, p^ + C^ 7^ the potential function vq as well as the 
integral wq diverge whilst these values remain finite for j > 0: 

vo = r^== / ^(a;'; ^i) dx + C(po), v, = C»(po), V2 = 0(po), 

PoVP^ + C -^-1 



Wo = , 3f? 

PoVp^+e 



Xi \x. 






with K being the Jacobian elliptic function K{m) = / 

Jo 



de 



Vl — "rn? sin 6 ' 



TTl 

2VP^ + C^ 



K'^'^^ is now on the other sheet of the Riemann surface. 
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Now, we define the curve Tu {X2 = —Xi): 
For all values Xl ^Tu, vq diverges as well: 






= 0(14) 



""^^^--W^'^''^'"^ w,thC(.Y.)^0 

In the limit po ~^ 0, this leads to finite solutions X^^^ = K^^^po (j = 1, 2) 
of Jacobi's inversion problem and eventually to 

f = ( p' r' ^^^^ .^,\ 1 

'''^Vv/p^+C^A(^) ^{x^-x(){x^-xi) ^^''J^ ■ 

2. The black hole limit: The above simultaneous limit is performed such that 

remains finite. We then get for po -^ (with p = rsin6, ( = rcosB): 

2Qu vricos^ _ vri 

r 2r^ 2r 

Again we can follow the procedure of pO[. Since Kj -^ {j = 1,2), the 
integrals in Jacobi's inversion problem (jl^) become elementary, and thus 
for / one obtains 

2rfi(7 — 1 — icos^ 
2rQu + 1 — i cos 6 ' 

i.e. the (r > part of the) extreme Kerr solution with J = 1/(4^2^) = M^. 
The constant Qu plays the role of the 'angular velocity of the horizon'. 

3. The non asymptotically flat ultrarelativistic lim,it: As shown in [^ for the 
rigidly rotating disks, our more general solutions also allow for a different, 
not asymptotically flat limit. To achieve this, one has to consider flnite 
values of r/po- A coordinate transformation 

(^ = ip-nut9 = 9, i= C{Xi) t (hence r/po = 2ffi[/) 



^^0 = — ^r:^— , vi = - — , V2 = 0. 



CiX,: 

yields a transformation to a corotating system (with angular velocity Qu) 
combined with a rescaling of r and t. In the limit po — > and C{Xi) —* 0, 
the resulting Ernst potential / [which is related to the Ernst potential /^ 
within the above corotating system by / = f[j/C^{Xi)] still describes a disk 
and is regular everywhere outside the disk. However, it is not asymptotically 
flat. 
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2.4 Backlund limit 

In the limit of real and negative values of Xf we obtain purely imaginary val- 
ues for Ki and K2 with Ki = K2. Then the reformulation (|l^) of the ex- 
pressions for / and Jacobi's inversion problem, leads to [with the abbreviation 

r, = ^{K^-iz)iK, + iz)] : 



f 




dK 



dK 



- [V2 - K^Vq] 



ITT 



KW {K-Ki)^{K-iz){K + iz) 



Vi + KiVq 

ri 



dK 



X(2) 



{K + Ki)^{K-iz){K + 



Vi - KiVq. 



\z\ 



-{v2-Kfvo) 



Applying the substitution X{K) = \J{K — iz)/{K + iz) 
system of equations for /, A'-^^A*^^\ and (A^^^ — A*^"^^) : 

awa(^) + (a(^)-aW)-i 

^ ~ A(i)A(2) - (A(2) - XW) - 1 - 

AWA(^)-Ai(A(^)-AW)-A^ 
AWA(2) + Ai(A(2)-AW)-A2 '' 
xmxm _ A*(A(2) - A(^) - Af 
AWA(2) + At(A(2)-A«)-Af 
[with A(^') = A(ir(^')),Ai = 
The solution for / is given by 



we get the following 



^r\{v\+K-i_vo) 



-^rl{vi-Kivo) 



A(i^i),At = l/Ai,rt = fi]. 



/ = /o 

where 
/o = e' 



1 1 1 

■I Aifti \\al 

1 \\ Af 



1 1 

Aitti X{a\ 

\l Af 



(15) 



-(v2-Klvo) 



ttl 



1 -exp[ri(t;i + KiVq)] 



ai 



1 



1 +exp[ri(fi + KiVq)]' 
Equation ( pIS]) represents a Backlund transformation of the real seed solution /o, 



see [|T0|, |TT1. As a consequence of (|^), ai satisfies the Riccati equations!^ 



«!,. 



OLi,z = Ai(ai - 1)77-^5 
fjo 

^°It is interesting to note that the special choice of the integration constant arising from (|§) 

guarantees the absence of a singular behaviour of the final solution at p = —^[Ki]. 



1/2 T\ fo,z 
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Hence, in this limit, our solutions are physically interesting Backlund transforms 
of nontrivial seed solutions of the Weyl class. 
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FIGURE CAPTIONS 

FIGURE 1: 

Outside the hatched region of the parameter space of Xf regular solutions have 
been found. Physical characteristics of solutions corresponding to values inside 
the allowed region are indicated. 

FIGURE 2: 

Examples for differentially rotating disks. The dimensionless quantities poCp, poQ 
and the function h are plotted against the normalized radial coordinate p/po and 
X, respectively, for 

(a) X2 ^ -1/2 + i/3 (here Xf G F,) (c) Xf = -3/2 + i/5 

(b) Xf = -2/3 + i/2 (d) Xf = -4. 



FIGURE 3: 

Contourplot of M'^/J. The stripes correspond to intervals of length 1/20; they 
are white for Newtonian disks {M'^/J -^ 0) and black for highly relativistic disks 
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Figure 1: Outside the hatched region of the parameter space of Xf regular 
solutions have been found. Physical characteristics of solutions corresponding 
to values inside the allowed region are indicated. 
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Figure 2: Examples for differentially rotating disks. The dimensionless quantities 
PoCpi Po^ and the function h are plotted against the normalized radial coordinate 
p/po and X, respectively, for 

(a) Xf ?« -1/2 + i/3 (here Xf e T^) (c) X^ = -3/2 + i/5 

(6) Xf = -2/3 + i/2 (d) Xf = -4. 
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Figure 3: Contourplot of M"^ /J. The stripes correspond to intervals of length 
1/20; they are white for Newtonian disks {M'^/J -^ 0) and black for highly 
relativistic disks {M'^/J -^ 1). 



